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Abstract 

The search for new physics in single- and multi-photon final states with large 
missing energy at LEP and future e + e~ colliders requires precise predictions for the 
Standard Model irreducible background. While at LEP1 the theoretical situation 
is under control, going to LEP2 (and beyond) some improvements are necessary. 
To approach the aimed 0(1%) theoretical accuracy, the tree-level matrix elements 
for the processes e + e~ — > vvnj, with n = 1,2,3, are exactly computed in the 
Standard Model, including the possibility of anomalous couplings for single-photon 
production. Due to the presence of observed photons in the final state, particular 
attention is paid to the treatment of higher-order QED corrections. Comparisons 
with existing calculations are shown and commented. An improved version of the 
event generator NUNUGPV is presented. 
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1 Introduction 



The production of one or more photons and missing energy in high-energy electron- 
positron (e + e~) collisions is a process of great interest for the scientific programme of 
LEP and future e + e~ linear colliders |TJ, 0. 

In the context of precision tests of the electroweak interactions on top of the Z res- 
onance (LEP1) P, ||, the Standard Model (SM) process e + e~ — > vvj, where v refers 
to light neutrinos and 7 to a detected, energetic photon, has been successfully used for 
the determination of the number of light neutrino species, in agreement with the result 
obtained via the measurement of the partial width of the Z-boson into invisible particles. 

Above the Z resonance, i.e. for current data taking in the LEP2 energy range (y/s ~ 
160-200 GeV) H] and for planned experiments at the TeV scale 0, the events with 
single- and multi-photon final states plus missing energy (i^) play an important role 
in the search for new phenomena beyond the SM [jl], 0. Actually, the SM processes 
e + e~ — > uPwy, with n = 1,2,..., are the largely dominating irreducible backgrounds 
to a New Physics (NP) signature consisting of one or more photon(s) and nothing else 
seen in the detector. Such events can indeed originate from various mechanisms, both in 
gravity- and gauge-mediated supersymmetric models as well as in scenarios with strong 
electroweak symmetry breaking || . In supersymmetric extensions of the SM, neutralinos, 
gravitinos and sneutrinos are the lightest supersymmetric particles (LSP) yielding the 
content of missing energy to the events. Another interesting example of new phenomenon 
giving rise to n.7 + 1$ final states is the production of a pair of fourth-generation, heavy 
neutrinos in association with initial-state radiation (ISR). Finally this signature can be 
useful to study anomalous triple and quartic gauge couplings. In particular, it can be 
conveniently used for the study of anomalous couplings as a complementary channel with 
respect to processes such as e + e~ — > 4/, 6/, without contamination coming from the 
vertices involving at least three massive gauge bosons. 

The present situation of the theoretical calculations for the above quoted processes 
can be considered as satisfactory for the purposes of data analysis at LEP1. Going to 
LEP2, the typical SM cross section is of the order of a few picobarn, yielding thousands 
of events collected in the four LEP experiments. Hence, there is a demand for theoretical 
predictions with an accuracy of the order of 1% for the rate of e + e~ — > vvwy events in 
the SM. Furthermore, as it will be discussed in the following, the situation concerning the 
theoretical calculations is not completely satisfactory. In particular, a careful treatment 
of higher-order QED corrections to processes with detected photons in the final state 
becomes mandatory for a meaningful comparison between data and theory. 

Given the above physics motivations, the aim of the present paper is to perform an 
exact tree-level calculation in the SM of the e + e~ — > vvw-f cross sections, with n = 1, 2, 3, 
supplemented with the most phenomenologically relevant and presently under control 
radiative corrections. While for one and two photons in the final state exact matrix 
element calculations have recently appeared in the literature, as it will be discussed in 
detail in the next Section, no results for three photons are available yet. Concerning 
ISR, particular care is devoted to the implementation of the effects due to (undetected) 
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photon emission. An improved version of the event generator NUNUGPV |§, based on 
the results here described, is presented for data analysis. 

The content of the paper is as follows. In Sect. 2 the status of the available calcula- 
tions and related programs is critically reviewed, putting special emphasis on the needs 
of improvement for current experiments at LEP. In Sects. 3 and 4 it is described how 
some of the existing approximations are overcome in the present study, providing various 
numerical results. The exact SM calculation of the lowest-order amplitudes for the sig- 
natures e + e~ — > vvny, with n = 1, 2, 3, is discussed in Sect. 3, while Sect. 4 is devoted 
to the implementation of higher-order QED radiative corrections. In passing, numerical 
results for integrated cross sections as well as comparisons with existing calculations are 
shown and commented. The impact of the theoretical improvements on exclusive photon 
distributions is discussed and compared with typical NP effects in Sect. 5. The main 
conclusions and open issues are drawn in Sect. 6. 

2 Status of the theoretical predictions and generators 

Concerning the process e + e~ — > described by the Feynman diagrams depicted in 
Fig. D, several calculations, with a different degree of accuracy, are known in the literature. 




Figure 1: The tree- level Feynman diagrams for the process e + e — > wy. 
Before the start of LEP/SLC operations, the status of the theoretical calculations and 
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related computational tools has been summarized in Ref. |T(| and, more recently, in view 
of the energy upgrade of the LEP collider from LEP1 to LEP2 regime, in Ref. [O . 

First attempts |L2| to compute the cross section associated to the graphs of Fig. [l| 
consider the limit My/ — > oo and neglect the contribution of the non-abelian WW'j vertex 
(the so-called Point Interaction Approximation). In other approximate calculations the 
invisible neutrino-pair cross section is dressed with some (universal) radiation factor to 
attach one external photon to the charged fermion legs, e.g. by using an angular dependent 
radiator fl3|, a parton shower (PS) algorithm (as in the program PYTHIA) |l5j or the 
Yennie-Frautschi-Suura (YFS) exclusive exponentiation (as done in KORALZ) |I6|| . By 
construction, these approaches allow to account for the leading contributions due to the 
collinear and infrared singularities but need to be corrected for the effect of sub-leading 
terms and/or internal photon radiation from the off-shell W boson (as recently done in 
an approximate way in Ref. [|T7|) that are contained in the exact matrix element. The 
first complete calculation of the matrix element of the process e + e~ — > vv^f was done in 
Ref. . The corresponding exact matrix element is implemented in the event generator 
MMM [ 19 1 . In Ref. [18|], by working in the approximation of neglecting in the squared 
matrix element terms with at least three boson propagators, also a compact, analytical 
expression for the differential spectrum on the energy and angle of the observed photon 
is obtained, yielding the result 

do 



a 



d cos "Q^dk 



12tt 



s'k 
sk+k„ 



(1) 



The meaning and explicit expression of the symbols entering eq. ([]]) can be found in 
Ref. [I8|. The photon spectrum of eq. ([]]) contains the bulk of the contributions due to 



iy-boson exchange and agrees within 1% with approximate calculations discussed above 
for center of mass (cm.) energies around the Z resonance |L4|]. The analytical photon 
spectrum of eq. (|1|) is implemented in the event generator NUNUGPV [[|. 

Concerning radiative corrections, the exact one-loop electroweak corrections to e + e~ — > 
vv^j process are not yet available. What is presently known is the full set of one-loop elec- 
troweak corrections to the "sub-process" e + e~ — > Zj, with an on-shell final Z boson |20 



and the subset of one-loop QED corrections to the Z-exchange contribution to z/z/7 final 
state [18 , [2~Tf . However, since, contrary to LEP1, the W^-boson contributions are essential 
at LEP2, as clearly shown in Fig. ^ for different selection criteria, the still missing part 
of the full one-loop calculation due to the corrections to the dominant diagrams with 
Vy-boson exchange is necessary in order to reach the aimed theoretical precision at the 
1% level. Actually, one can see from Fig. [|that the relative contribution of iy-diagrams 
(squared modulus and interference terms) can amount to about 60% (80%) of the full 
cross section, including or excluding, via a cut on the missing mass of the event, the Z 
return. Therefore, in the absence of a full 0(a) electroweak calculation, the goal of a 1% 
theoretical accuracy in the predictions for single-photon plus missing energy production 
turns out to be difficultly reachable at present. However, in order to take care of the most 
sizeable higher-order corrections, the lowest-order calculations are typically improved by 
the inclusion of the (large) effects due to ISR. In mostly used computational tools such a 
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Figure 2: The relative size of the ^-exchange diagrams with respect to the total cross 
section of the signature for different experimental cuts. The quantity Ao"i 7 is defined 
as: A<7i 7 = 1 — 3<j(t' Ai z/ A1 7)/cr(z/z/7). In Fig. 2a the four lines correspond to the cuts on 
the photon energy and angle specified in the plot; Fig. 2b is the same as Fig. 2a, with an 
additional cut on the missing mass of the event. 



contribution is taken into account via traditional algorithms for computing QED radiative 



corrections in the leading logarithmic (LL) approximation, such as the PS algorithm |22 



(as in PYTHIA), the Structure Function (SF) approach |23| (as in MMM and NUNUGPV, 
but also in PYTHIA) and YFS exclusive exponentiation |24| (as in KORALZ). Two dif- 
ferent variations of the SF method are employed in MMM and NUNUGPV. SFs in strictly 
collinear approximation are used in MMM, while in NUNUGPV pt-dependent SFs |13|. [25[ 
are implemented to improve the treatment of ISR with the effect of the transverse degrees 
of freedom at the LL level || (this point will be further discussed later). 

The previously quoted programs KORALZ, MMM and NUNUGPV are the standard 
Monte Carlo generators used by the LEP collaborations for the analysis of the data relative 
to the events e + e~ — * 1/1/7(7) JJ. As already emphasized, these programs differ both in the 
treatment of the lowest-order matrix element and in the implementation of higher-order 
corrections. Actually, the agreement between the above generators, whenever compared 
to estimate the theoretical precision against the statical accuracy in the measured cross 
section (of the order of a few per cent), is fairly good for exclusive single-photon final 
states (say for photon energies above approximatively 10 GeV), whereas for low photon 
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energies and more inclusive final states the agreement is not at present satisfactory |26 



Therefore, since in the search for new stable, neutral particles one is interested in detecting 
an excess of events at low photon energies [p6fl , the present status points out the need of 
improving theoretical predictions to avoid a loss of sensitivity to NP searches in radiative 
events at LEP2. 

Let us come now to discuss the present status of the theoretical predictions for the 
process e + e~ — > z/z/77, which is the most relevant SM background to a signature with 
two acoplanar photons and large missing energy. Only very recently, dedicated calcula- 
tions appeared in the literature, as motivated by the search of anomalous 77 + Ifi events 
at LEP. A complete diagrammatic calculation, using the helicity amplitude technique, 
supplemented with collinear SFs to account for ISR, was done in Ref. [PTfl . This paper 
confirmed a previous evaluation in Born approximation contained in Ref. |p8 |. Approx- 
imate predictions for the process of interest are obtained by the LEP collaborations by 
using the above quoted programs with QED "dressing" of the neutrino-pair cross section, 
namely PYTHIA (via the PS) and KORALZ (via the YFS method). Further, modern 
packages for the automatic calculation of Feynman amplitudes, such as GRACE [29| and 
CompHEP |3"0 ], are used by the experiments to calculate the cross section and gener- 
ate events. Both packages implements collinear SFs for ISR, with an option for PS in 
GRACE. 

Quite recently, an extensive comparison, at the level of total and differential cross 
sections, between all the available calculations has been performed in Ref. 



311. This 



detailed comparison shows a (dis) agreement between independent calculations at 10-20% 
level (or worse) and, more generally, an unsatisfactory situation about the software for 
the analysis of the data with acoplanar photons (see also last paper in Ref. Pf). 



3 Calculation of the lowest-order cross sections 

In order to approach the 0(1%) theoretical precision and improve the predictions of the 
earlier version of the program NUNUGPV, the lowest-order matrix elements associated 
to the processes e + e _ — > vvwy (with n = 1,2, 3) have been exactly calculated in the SM. 

The matrix element for single-photon production has been computed by means of 
helicity amplitude techniques P2f , including the possibility of anomalous Afc 7 and A 7 
contributions to the WW'j coupling. This offers the possibility of exploiting the LEP 
statistics relative to 7 + 1$ events in order to put constraints on such anomalous terms. 
The diagrammatic calculation has been cross-checked by using the algorithm ALPHA |33|] 
and found to be in perfect agreement. The three-body phase space has been generated 
recursively, by decomposing the phase-space element c?$ 3 as follows 

d$ 3 (P -> k, q h q 2 ) = rf$ 2 (P -> k, Q) (2ir) 3 dQ 2 d<S> 2 (Q -> q h q 2 ), (2) 

where P = pi + p 2 is the total incoming momentum, Q = q\ + q 2 is the total momentum 
carried by the neutrinos and k is the photon momentum. The independent kinematical 
variables are chosen to be: E y , cos 1 ^, 7 , cos^*, (f>*, where the photon variables are 
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Figure 3: The relative deviation between the exact e + e~ — > vv^ cross section and the ap- 
proximated one, as obtained via eq. (1) implemented in the earlier version of NUNUGPV. 
The differences are shown for the event selections given in Fig. 3a, with an additional cut 
on the missing mass in Fig. 3b. 



generated in the cm. frame, while the neutrino ones in the reference frame where ql + (j2 = 
0. Since the photon energy can be expressed in terms of the invariant mass Q 2 via the 
linear relation E y = (s — Q 2 )/2y / i, where s is the total cm. energy, it is convenient to 
generate Q 2 in such a way to sample directly the leading matrix element configurations, 
that are due to the emission of a soft photon or a hard, Z return one. This importance 
sampling strategy is followed in the Monte Carlo integration in order to cure the variance of 
the matrix element in correspondence of the infrared and Z return peaking behaviour. The 
photon angle is generated according to the weight function p(cos $ 7 ) oc 1/(1 — /? 2 cos 2 $ 7 ) , 



with (3 = J 1 — 4m 2 /s, to take care of the collinear peaking. 

The exact treatment of the single-photon matrix element upgrades the released version 
of NUNUGPV, based on the photon spectrum of eq. (ffl) of Ref. |T8L to include previously 
neglected jy-boson effects relative to contributions with at least three boson propagators. 
The size of such previously neglected effects is shown in Fig. [|, for typical event selections 



used by the LEP experiments [g, ^6[. The calculation of the single-photon cross section 
obtained with the exact matrix element is compared with the cross section resulting from 
the integration of the photon spectrum of eq. ([!]). The relative difference between the two 
calculations is at a few per cent level, both without and with a cut on the missing mass, 
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in agreement with the degree of approximation stated in Ref. [ 18fl . However, it should be 
noticed that an exact treatment of the lowest-order matrix element is actually mandatory 
at LEP2 if a theoretical accuracy of the order of 1% is aimed at. 
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Figure 4: The tree-level cross section for the process e + e~ — > z/z/7 as compared with 
the cross section with higher-order QED corrections, obtained by using collinear SFs. 
Two typical selection criteria, specified in Fig. 4a, are considered, including (Fig. 4a) and 
excluding, via a cut on the missing mass (Fig. 4b), the Z radiative return. 



The already quoted algorithm ALPHA, that is conceived for the automatic compu- 
tation of tree-level multi-particle production amplitudes without any need of Feynman 
graphs expansion, has been employed for the calculation of the matrix elements with two 
and three photons in the final state. For the case zyz/777 the calculation here presented 
is the first one appearing in the literature. Without entering the details of the algorithm 
ALPHA, which is fully discussed in the literature [Q, it is worth noticing, for the aim 
of the present study, that the predictions of this automatic algorithm have been already 
compared with the diagrammatic results for the processes e + e~ — *> 4/ [IT, 134], showing ex- 



cellent agreement, and also successfully used to obtain original results for other reactions 
such as 77 — > 4/ j35|, confirmed in Ref. p6[ , 



e 1 e 



4/ + 7 PH an d e + e 



,6/ g. 

Concerning the treatment of phase space, the decomposition introduced for two photons 
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in the final state reads as follows 

rf $ 4 (P ^ kl , k 2 , qx, q 2 ) = d$ 3 (P -> h, k 2 , Q)(2n fdQ 2 d^ 2 (Q -> g 1; g 2 ), (3) 

where, as before, Q = qi +q 2 is the total momentum carried by the neutrinos and ki, k 2 are 
the momenta of the two photons. The independent kinematical variables are chosen to be: 
Q 2 , E y> i, E 1>2 , cos^ 7i i, 7i i, 7i i2, cos6**, 4>*, where the photon variables are generated 
in the cm. frame, while the neutrino ones in the reference frame where q[ + q~2 = 0. 
As in the case of single-photon production, the energies of the photons are generated 
according to the soft and Z return peaking structure, whereas one of the photon angles 
follow the collinear behaviour. A generalization of the above strategy is also followed for 
three photons in the final state. 

Some numerical results for the lowest-order cross sections of the processes e + e" — > vvj 
and e + e~ — > z/z/77 are shown, as functions of the cm. energy in the LEP2 energy range, 
in Fig. f| and Fig. |. The input parameters used throughout the present study are the 
same as in Ref. j!4l|. By comparing Fig. Eland Fig. & it can be seen that the cross section 
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Figure 5: The same as Fig. 4 for the process e + e — > z/z/77. 

for the signature ^z/77 is about a factor 10-100 smaller than the cross section for z/z/7, 
the reduction factor being strongly dependent, as expected, on the imposed photon cuts. 
A similar ratio is present in the LEP2 range for the z/z/777 cross section with respect to 
the one for the ^z/77 final state. 
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4 Treatment of higher-order QED corrections 



Because the level of accuracy demanded to the calculations for single- and multi-photon 
production with missing energy should reach the 0(1%) precision, the most important 
radiative corrections must be necessarily included. To this end, all the calculations known 
in the literature take into account the effect of ISR. As already discussed in Sect. 2, 
this goal is achieved by using standard algorithms for universal photonic corrections, 
such as the PS method, the YFS exponentiation and the SF approach. In particular, 
the SF approach, because of its simplicity, is certainly the most widely used algorithm, 
implemented in many generators of interest here, such as CompHEP, GRACE, MMM 
and NUNUGPV. More precisely, as discussed in Sect. 2, all the programs make use of 
SFs in strictly collinear approximation, while in NUNUGPV ^-dependent SFs [|T^, 
are implemented to improve the treatment of ISR by including pt/pL effects. Because of 
the presence of photons among the observed final-state products, the inclusion of the ISR 
requires a particular care. This caution is further motivated by the very large enhancement 
of the lowest-order cross sections as due to the ISR and clearly visible in Fig. £| and 
Fig. |5|. As can be seen, the single- and double-photon cross section are significantly 
enhanced when considering typical event selections. The enhancement factor is about 1.3 
when including the Z return and about 2 when excluding it. As done in many practical 
applications [14], [1^, p7|, the above results can be simply obtained by convoluting the 



hard-scattering cross section of interest with collinear SFs, according to the factorized 
formula 

Ccoii = J dx\dx2 D(xi, s)D(x2, s) daQ(cuts). (4) 

It allows to take into account the impact of higher-order QED corrections, due to photon 
emission before the hard-scattering reaction (pre-emission), at the LL level. It has been 
checked that the dominant effects of the ISR obtained according to eq. (§) numerically 
agree, for different experimental set up, with those already known in the literature for 



one and two photons in the final state [14, 27, Equation (m is a good approximation 



to QED radiative corrections in the LEP1 energy regime at the LL level, since, with the 
standard selection criteria, hard "pre-emission" photons imply the production of vv pairs 
off Z-resonance and are therefore inhibited. On the contrary, going to LEP2 it can be 
easily realized that the implementation of the ISR as given by eq. (f|) is an approximation 
that clearly fails whenever the photonic degrees of freedom, of the pre-emission and hard- 
scattering process, respectively, overlap in the same phase space region. Actually, because 
the collinear SFs can be seen as the result of an integration over the angular variables 
of the photon radiation, eq. (^) does not take into account the correct statistical factor 
to be included for the presence of identical particles in the final state. Furthermore, if 
the pre-emission photon is detectable, the reconstruction of the event via eq. (f|) is only 
approximate and this might imply an additional inaccuracy. Therefore, one should expect 
that the implementation of the ISR as given by eq. (^) leads to an overestimate of the 
higher-order QED corrections. This effect is clearly dependent on the photon(s) detection 
criteria and can be expected to be not negligible with respect to an 0(1%) theoretical 
accuracy. An estimate of the effects due to the phase space overlapping of the IS pre- 
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emission photons with the observed ones can be obtained by supplying the QED SFs with 
the transverse degrees of freedom. Actually, the generation of the angular variables at the 
level of the ISR gives the possibility of rejecting in the event sample those pre-emission 
photons above the minimum detection angle, thus avoiding "overlapping effects". 

According to such a procedure, the cross section with higher-order QED corrections 
can be calculated as follows (for the realistic data sample of at least one photon) 

(j l7(7) = J dx^dc^dc^ D(x u c { ^;s)D(x 2} c {2) ;s)Q(cuts) 

x (da 1 ^ 1 + da 2 ^ 1 + da^ + ..) , (5) 

where D(x, c 7 ; s) || is a proper combination of the collinear SF D(x, s) with an angular 
factor inspired by the leading behaviour l/(p-k). The latter is introduced to generate the 
angular variables of the pre-emission photons. According to eq. d^), an "equivalent" pho- 
ton is generated for each colliding lepton and accepted as a higher-order ISR contribution 
if: 

• the energy of the equivalent photon is below the threshold for the observed photon 
E^^min, for arbitrary angles; or 

• the angle of the equivalent photon is outside the angular acceptance for the observed 
photons, for arbitrary energies. 

Within the angular acceptance of the seen photon(s), the cross section is evaluated 
by summing the exact matrix elements for the processes e + e~ — ► vvwy, n = 1,2,3 
{da 1 " 1 , da 21 , da 3 " 1 ) . Notice that from the point of view of computing a 17 '' 7 -' the real con- 
tributions da 111 , n > 2, represent the "hard" radiative corrections to be matched with 
the universal soft + virtual ones accounted for by the SFs. Therefore they are in prin- 
ciple necessary at all orders. The truncation of hard radiative corrections at the level 
of da 37 introduces a spurious infra-red sensitivity in radiative corrections at the order 
a 4 In (E I Ej !m i n ) which, from the practical point of view, is completely negligible at real- 
istic E lim i n . Since the radiative corrections implemented by means of this procedure are 
at the LL level, its theoretical error is dominated by missing truly 0(a) corrections. 

In order to quantify the overestimate introduced by the collinear SFs in the calculation 
of the ISR via eq. @, the relative difference between eq. (f|) and eq. (H) is shown in Fig. ||, 
for several photon detection criteria. As can be seen, at LEP1 the overlapping effects are 
contained within a few per mille, and therefore negligible on the scale of the experimental 
accuracy. Going to LEP2 energies, the impact of "overlapping effects" [] varies within 
1-4% when including the Z return and is still larger, reaching 10%, when imposing a cut 
on the missing mass, as usually done in realistic event selections. This effect is therefore 
important in the light of the aimed theoretical precision. It should be noticed, however, 

1 Notice that, as already remarked, by "overlapping effects" two effects are understood, which do 
occur simultaneously whenever two or more photons are detected: in eq. (|J) i) the statystical factor is 
incorrectly accounted for, ii) the reconstruction of the event is approximate. 
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Figure 6: Contribution of "overlapping effects" (see the text for definition) to the cross 
section for the process e + e~ — > 1/1/7(7). In Fig. 2a the four lines correspond to the cuts 
on the photon energy and angle specified in the plot; Fig. 2b is the same as Fig. 2a, with 
an additional cut on the missing mass of the event. 

that, where the difference reaches a ten per cent size, the overall effect of ISR is to 
enhance the tree-level cross section of a factor of two. A qualitative explanation of the 
overlapping effects discussed above can be given as follows. The overestimate of radiative 
corrections takes place when the pre-emission photon can reach the observability region 
for the detected photon. At LEP1, and with standard selection criteria, the emission of 
multiple detectable photons implies that neutrino production occurs off zT-resonance and 
therefore the overlapping effect is naturally suppressed by the dynamics. At LEP2, where 
this suppression is no longer active, the overlapping effects can become more sizeable, 
depending on the angular acceptance and minimum energy of the observed photons, and, 
more generally, in the presence of additional cuts on the four-momenta of the observed 
photons (as in the case of a cut on the missing mass). 

Analogously to eq. (0), the QED corrected cross section for the signature of at least 
two-photons in the final state can be cast as follows 

a 27(7) = J dx x dx 2 dc { ^ ) dc^ ) D(x l ,c^;s)D{x 2 ,cf ) ;s)Q{cuts) 

x (da 2 ^ + da 3 ^ + ...) . (6) 
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It is worth noticing that da 31 in eq. (H) plays the same role as da 21 in eq. @, and hence is a 
key ingredient when considering the signature with at least two photons in the final state. 
Numerical results for such a signature are given in Fig. [j], showing the relative difference 
between eq. (|j) and eq. Also for the process e + e~ — > 2/^77(7) the implementation of 
the ISR via collinear SFs can lead to an overestimate of the physical cross section of the 
order of several per cent. 
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Figure 7: The same as Fig. 6 for the process e + e — > 2/^77(7). 



5 Distributions 

The results shown in the previous sections refer to integrated cross sections. However, 
in the search for NP effects, one is also interested in more exclusive distributions involv- 
ing the photon(s) energy and angle. Therefore, it is worth studying how the theoretical 
improvements discussed above on tree-level matrix elements and higher-order QED cor- 
rections have an impact on photon distributions with respect to typical NP deviations. 

To this aim a few results obtained by using the improved version of NUNUGPV as 
an event generator are illustrated in Figs. 3-|l0|. Figure |8] and Fig. |] refer to the signature 
e + e~ — > ^7(7), while Fig. |9] is relative to the process with at least two photons, i.e. 
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Figure 8: "Soft" p] distribution for the leading photon of the process e + e~ — > 1/1/7(7) at 
yfs = 190 GeV. The cuts on the observed photons are E 7 > 1 GeV and 12° < $ 7 < 168°. 
Exact tree-level calculation (solid line) versus: approximate analytic spectrum of Ref. flTf] 
(dashed line) (Fig. 8a), exact SM calculation plus anomalous couplings contribution 
(dashed line) (Fig. 8b), exact SM calculation plus massive (m„ = 50 GeV) fourth family 
neutrinos (dashed line) (Fig. 8c). In Fig. 8d ISR correction is implemented according to: 
eq. (H) (solid line), eq. ([|) (dashed line) and via SFs with pt/VL effects as in the released 
version of NUNUGPV (dotted line). 



e + e~ — > 1/1/77(7). For the sake of comparison, the histograms shown in Figs. |q-|10| are 
normalized to the same luminosity. A typical LEP2 energy of ^fs = 190 GeV is considered, 
with the cuts E 1 > 1 GeV and 12° < $ 7 < 168° for the observed photons. 

As an example, the p] distribution of the most energetic photon (leading photon) is 
shown in Fig. |8| for the "soft" region and in Fig. |9] for the "hard" one. The first three 
plots (Figs. 8a-8c, Figs. 9a-9c) do not take into account the contribution of ISR (Born 
approximation), while the last one (Fig. 8d, Fig. 9d) shows QED corrected distributions 
according to different realizations. The aim is to show how typical NP effects may compete 
with an improved calculation of the SM background, both at the level of tree-level matrix 
elements and higher-order QED corrections, previously discussed. Actually, in the first 
three plots the solid line is the lowest-order prediction obtained by means of the exact 
single-photon matrix element as compared with the following results (dashed histograms) 
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Figure 9: The same as Fig. 8, in the "hard" p[ region. 

1. approximate lowest-order photon spectrum of eq. (|I|), implemented in the released 
version of NUNUGPV (Fig. 8a, Fig. 9a); 

2. exact single-photon matrix element with anomalous WW'y coupling, corresponding 
to the parameters choice A/c 7 = A 7 = 5 (Fig. 8b, Fig. 9b); 

3. exact single-photon matrix element with additional contribution due to the produc- 
tion of a pair of fourth-generation neutrinos, with standard couplings and mass of 
50 GeV (Fig. 8c, Fig. 9c). 

As can be seen, an exact treatment of the lowest-order matrix element is mandatory in 
order to obtain fully reliable exclusive photon distributions and to avoid loss of sensitivity 
in the search of (small) NP deviations. Also a careful formulation of ISR is necessarily 
required for an appropriate simulation of photon distributions, as shown in Fig. [8]d and 
Fig. |]d. In these plots the three following different implementations of higher-order QED 
corrections are compared: improved treatment of ISR via eq. (||) (solid histogram), sim- 
ulation of ISR via collinear SFs as given by eq. ([3]) (dashed histogram), treatment of ISR 
via SFs with pt/pi effects as in the released version of NUNUGPV (dotted histogram). 
The differences between the three implementations of ISR are certainly comparable or 
larger than the deviations introduced by NP, clearly illustrating the need of a proper 
treatment of pre-emission QED effects. In particular, the collinear SFs lead, as expected, 
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to an overestimate of events in the soft-photon region, while the implementation of ISR 
via SFs with pt/pL effects, as in the released version of NUNUGPV, is responsible for an 
excess of events in the large p] region, as a consequence of the approximate treatment of 
p t contributions outside the collinear approximation. 




Figure 10: E 1 and pi distributions for the leading (Figs. lOa-lOb) and next-to-leading 
photon (Figs. lOc-lOd) of the process e + e~ — > 1/^77(7) at yfs = 190 GeV. The cuts 
on the observed photons are E 1 > 1 GeV and 12° < $ 7 < 168°. ISR correction is 
implemented according to: eq. (|]) (solid line) and eq. ([|) (dashed line). 

A comparison between photon distributions obtained by means of different simulations 
of ISR is also shown in Fig. |H] for the e + e~ — > 1/1/77(7) events. The implementation of 
ISR via collinear SFs as in eq. (|]) (dotted histograms) is compared with the improved 
formulation of eq. (|6]) (solid histograms). The energy and p t distributions of the leading 
(Figs. lOa-lOb) and of the next-to-leading (Figs. lOc-lOd) photon are considered. It can 
be seen that the implementation of ISR via collinear SFs, as done in practice in many 
computational tools for 1/1/77 final state, lead to an excess of events in the soft-photon 
region, as already noticed for the uP^ signature. Since the soft-photon region is of primary 
interest for the search of NP effects, a careful treatment of ISR, as given by eqs. (|§)-(§P, 
is essential to obtain actually precise predictions for the SM irreducible background. 
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6 Concluding remarks 



The search for new physics in single- and multi-photon final states with large missing en- 
ergy at LEP and future e + e~ colliders requires the best knowledge of the SM irreducible 
background. To this aim, precise calculations (and related computational tools) for the 
processes e + e~ — > uuny are presently demanded. Towards such a direction, the tree-level 
matrix elements for the SM processes with neutrino pairs and up to three photons in the 
final state have been calculated without any approximation. The exact treatment of the 
lowest-order transition amplitudes has been seen to be actually necessary in view of an 
expected precision at the 1% level. At this accuracy level, also a careful treatment of the 
(large) effect of the higher-order corrections introduced by ISR is unavoidable. Indeed, it 
has been shown that the usual implementation of ISR via collinear SFs, which is a good 
approximation at LEP1 energies and with the usual selection criteria, can lead to an over- 
estimate of the physical cross section much larger than 1%. In particular, in the presence 
of a missing mass cut, the ISR overestimate of the integrated cross section reaches the 
ten percent size. Furthermore, significant effects, whenever compared with typical NP 
deviations, have been shown to be present also in the photon distributions and should 
therefore carefully considered in a sensible experimental analysis. The remaining uncer- 
tainty in the present study is left to the yet unknown exact 0(a) electroweak corrections 
to the process e + e~ — * uuj. Such a complete calculation should be actually desirable to 
reach a theoretical error not exceeding the 1% level. 

As a result of the present study, an improved version of the event generator NUNUGPV 
is by now available. It includes the exact SM matrix elements for vvn^ production, with 
n — 1, 2, 3, and a careful treatment of the ISR including p t /PL effects. For single-photon 
production the possibility of studying the effects of anomalous couplings is included. 
Predictions for the production of a pair of hypothetic massive neutrinos can be also 
obtained. A sample of numerical results showing the potentials of the new version of 
NUNUGPV has been shown, with particular emphasis on the impact of the discussed 
improvements on the integrated cross sections and more exclusive distributions. The 
program can be used for a full analysis of single- and multi-photon events with missing 
energy at LEP2 and future high-energy e + e~ colliders. 
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